Temperature dependence of the upper critical field of an anisotropic 
singlet superconductivity in a square lattice tight-binding model 

in parallel magnetic fields 
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Upper critical field parallel to the conducting layer is studied in anisotropic type-II superconductors 
on square lattices. We assume enough separation of the adjacent layers, for which the orbital pair- 
breaking effect is suppressed for exactly aligned parallel magnetic field. In particular, we examine 
the temperature dependence of the critical field H C (T) of the superconductivity including the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO or LOFF) state, in which the Cooper pairs have non-zero center- 
of-mass momentum q. In the system with the cylindrically symmetric Fermi-surface, it is known 
that H C (T) of the d-wave FFLO state exhibits a kink at a low temperature due to a change of the 
direction of q in contrast to observations in organic superconductors. It is shown that the kink 
disappears when the Fermi-surface is anisotropic to some extent, since the direction of q is locked 
in an optimum direction independent of the temperature. 



Recently, upturn of the curve of the upper critical field 
as a function of the temperature has been observed in 
organic superconductors when the magnetic field is ap- 
plied in the directions parallel to the conducting lay- 
ers . This behavior has often been discussed [fz| — jl3|] in 
connection with a possibility of the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state p4Ul3]. It is known by many 
theoretical calculations |T^J17| ,|7j-p[ that the parallel crit- 
ical field exhibits upturn at low temperatures for the 
FFLO state in some models. However, in d-wave pairing, 
the curve of the critical field has a kink above the tem- 
peratures at which the upturn occurs when the orbital 
pair-breaking effect is negligible ]l8| , pd] , ^9|p3| , althought 
the kink has not been observed yet. In this paper, we 
show that the kink is suppressed by anisotropy of the 
Fermi-surfaces of some extent. 

Firstly, we briefly review theories and experiments 
which are concerned with the present theory. The su- 
perconductivity in type II superconductor is suppressed 
by applied magnetic fields due to both the Lorentz force 
and Zeeman energy for the singlet pairing. However, the 
Lorentz force is not effective, when the magnetic field 
is applied in directions parallel to the conducting layers, 
since the spatial motion of electron is restricted in the 
conducting layers. In this situation, the pair-breaking is 
mainly due to the Zeeman energy, and the upper critical 
field is known as the Chandrasekhar and Clogston limit 
(Pauli paramagnetic limit) pt^pjj . 

For example, for s-wave pairing, the critical field at 

zero temperature Hp^ is estimated from the zero field 

transition temperature Tc°^ by a formula [loHp^ « 1.25 x 

&bTc°^ in the weak coupling limit, where is the elec- 
tron magnetic moment. The critical field is maximum 
at T = with [d# P (T)/dT] T=0 = 0, and it is upwards 
convex at low temperature (d 2 Hp(T)/dT 2 < 0). 

Fuldc and Fcrrell and independently Larkin and 
Ovchinnikov |L5fl predicted that the critical field can ex- 



ceed the original Pauli paramagnetic limit by a pairing 
state with a non-zero center-of-mass momentum q. This 
state is called the FFLO or LOFF state, and as a di- 
rect consequence of the non-zero center-of-mass momen- 
tum, the superconducting order parameter oscillates in 
space. Larkin and Ovchinnikov found that the order pa- 
rameter has the spatial structure like cos(q • r) near the 
critical field in three dimensions rather than the struc- 
ture like exp(iq • r) by calculations of the free energies of 
the states, although both structures have the same upper 
critical field within the second order transition. 

The critical field of the FFLO state in three dimen- 
sions has a similar temperature dependence to the Pauli 
limit as mentioned above, wherease the value is slightly 
larger than the Pauli limit. Furthermore, the orbital 
pair-breaking effect strongly suppresses the FFLO state. 
The condition for the coexistence of the FFLO state and 
the vortex state was examined by Gruenberg and Gun- 
ther |2^] . When the FFLO state coexists with the vortex 
state, the spatial oscillation of the order parameter due 
to the FFLO state is along the vortex lines. 

On the other hand, in the two dimensions, the critical 
field as a function of the temperature is downwards con- 
vex (d 2 H c2 {T)/dT 2 > 0) at low temperatures 
This behavior, which is called upturn, is due to a Fermi- 
surface effect analogous to nesting effects of spin and 
charge density waves 0,^,^,|3|. In the two dimensions, 
the Fermi-surfaces of up and down spin electrons touch 
on lines by a shift by q at T = 0. It was found that by 
a calculation of the free energy that the order parameter 
has various two dimensional structures at low temper- 
atures and high fields in the two dimensions p^| . The 
origin of the two dimensional structure is the high crit- 
ical field, for which a gain in the polarization energy is 
more important than a loss in the condensation energy. 

Temperature dependences similar to the theoretical 
prediction including the upturn have been observed in 
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the organic superconductors |,@,§-§. Hence, the FFLO 
state has often been discussed as a candidate to ex- 
plain the critical field, although dimensional crossover 
effects in triplet and singlet superconductors are also 
possible mechanisms of the upturn |25|]2l| . How- 
ever, the observation of the lower critical field in k- 
(BEDT-TTF) 2 Cu(NCS) 2 §, a change in the temper- 
ature dependence of upper critical field due to the 
purity of the samples of A-(BETS) 2 GaCl 4 @, and 
that due to tilt angle of the magnetic field in k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br seem to support 

the FFLO state in these organic superconductors. In 
the phase diagram of A-(BETS)2FeCi4, the lower crit- 
ical field of the field induced superconductivity shows 
upturnlike behavior (d 2 H c2 (T)/dT 2 < 0) which is nat- 
urally explained by the FFLO state in combination with 
the Jaccarino and Peter mechanism |2^-|2{J. 

However, when we consider the FFLO state in the 
organic superconductors and the high-T c cuprates, we 
should note that d-wave pairing is an important can- 
didate of the pairing. The upturn of the critical field 
occurs also for the d-wave pairing as studied by Maki 
and Won and many other authors JTT| , p^|Jl^ ] . How- 
ever, in addition to the upturn at low temperatures 
(T < 0.06T C (0) ), the upper critical field H C (T) exhibits 

a kink at T w 0.06Ti°' ) due to a change in the direc- 
tion of q ]l8| , pl] , ^9| , p^[ . When we assume d x 2_ y 2-wave 
pairing, q is oriented to one of the crystal axes {J/x- 
axis or/y-axis) for lower temperatures, whereas it is ori- 
ented to the direction of (1, 1, 0) for higher temperatures 

(0.06T C (0) < T < 0.56T c (O) ). This behavior does not seem 
to be consistent with the experimental data of the upper 
critical fields §0,|-§]. 

Here, we should add a discussion for the case in which 
the vortex state coexists with the FFLO state in the pres- 
ence of the orbital pair-breaking effect. In this case, q of 
the FFLO state is fixed in the direction of the vortex line, 
and hence the kink does not occur. However, the tem- 
perature dependence of the critical field would change 
completely depending on the direction of the magnetic 
field within the directions parallel to the layers. 

In this paper, we examine the case in which the orbital 
pair-breaking effect is negligible and q is freely oriented 
to the optimum direction at each temperature and field. 
The orbital pair-breaking effect is negligible, when sep- 
aration of adjacent layers is sufficiently large [p5[ , the 
layers are sufficiently thin [^0|, and the magnetic field is 
exactly aligned in a direction parallel to the layers . 
It is shown that anisotropy of the Fermi-surface locks 
the optimum direction of q, and thus the kink of H C (T) 
disappears. We examine a tight-binding model on the 
square lattice ]g|,|26|] as an example of an anisotropic sys- 
tem. By change of the hole concentration, the shape of 
the Fermi-surface can be controlled between square and 
circular shapes. 

At zero temperature, the hole concentration depen- 
dence of the critical field has been examined in the square 



lattice tight-binding model ||. At hole concentration 
per a site ss 0.630, (i.e., electron density per a site 
n ~ 0.370), the nesting is most effective, and the critical 
field is extraordinarily enhanced near this hole concen- 
tration. The Pauli paramagnetic limit Hp was also cal- 
culated in this paper, and it was found that the FFLO 
critical field exceeds Hp" 1 at all hole concentrations ex- 
cept in a small region near the half-filling. 

In the present theory, we concentrate ourselves on the 
temperature dependence of the upper critical field. We 
consider sufficiently weak coupling of the pairing interac- 
tion and assume implicite small intcr-layer electron hop- 
ping, which makes the mean field treatment appropriate. 

The gap equation that we examine is 



(i) 



where /(e) is the Fermi distribution function, 7 Q (k) 
is a symmetry factor of the gap function, E^ = 
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/io|H| J?j, [io||n|]13[ |. By an approximation of the standard 
weak coupling theory, we obtain 
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and vp is the Fermi velocity. Here, p(e, 9) is the angle 
dependent density of states which satisfies 



k 



(4) 



for any function f (ek), and kp(9) is the magnitude of 
the Fermi momentum in the direction of 9, which is the 
angle measured from the fc^-axis. By solving equation 
eq. (2) with q fixed, we obtain a transition temperature 
T(h, q) or a critical field h(T, q) for the value of q. The 
final result of the transition tempreature is given by 



T c (h) =max(T(/i,q)) 



(5) 



As we mentioned above, we use the square lattice tight- 
binding model 



fk 



-2t (cos k x + cos k y ) — [i 



(6) 
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to express various anisotropic Fermi-surfaces. The cylin- 
drically symmetric Fermi-surface is obtained in the limit 
of [i — > — At, whereas a square Fermi-surface is obtained 
in the limit of fi — » 0. Figure 1 shows the shape of the 
Fermi-surfaces of fj, = — 0.5t, — 1.8t and — St. In this 
model, £ is written as C = ^ (qx(9, <pq) — 1) with q = tq/h 
and x(9,ifq) = [cos <,o q sin(fc x ) + sin^q sin(fcj,)]|k|=fc p (0), 
where <^ q is the angle between q and fc x -axis. We exam- 



ine d. 



2_„2-wave pairing 



7d 



, (k) = [cosfc^ - cosfc a ]| k | =fcF( 



(7) 



which is more favorable than d xy -wave pairing, because 
of the peak in the angle dependent density of state p(e, 9) 
near the saddle points of the electron dispersion relation. 
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FIG. 1. The solid, dotted and dashed lines show the 
Fermi-surfaces at fi — — 0.5t, — 1.8t and 3t, respectively. The 
lattice constant is taken as unity. 

Figure || shows the critical field for // = — 2>t, i.e., 
n ks 0.17. In the vertical axis, h = no\H\ is scaled by 
Ado = 2wd exp(— l/VNd{0)), which is rather different 
from the zero field BCS gap for anisotropic pairing 
As shown in Fig. |l|, the Fermi-surface has a nearly cylin- 
drical shape at this value of fi. The behavior of the criti- 
cal field is similar to that for the cylindrically symmetric 
Fermi-surface, which corresponds to p, — * —At. The op- 
timum direction of q of the FFLO state changes at a 
temperature T w O.lTi ', where the kink appears. Weak 
anisotropy does not change this behavior. It is confirmed 
by numerical calculations that the other directions q than 
ifq = and 7r/4 give lower critical fields than the highest 
one. 



FIG. 2. The solid and dotted lines show the critical fields 
for ipq — and n/4, respectively, when fi = —St, i.e., n « 0.17. 
In the inset, the solid and dotted lines show the tempera- 
ture dependences of q — |q| along the transition lines, for 
ifq = and 7r/4, respectively. The dashed and dot-dashed 
lines show the results of cylindrical symmetric Fermi-surface, 
(i.e., — > —it), for ipq = and 7r/4, respectively. 

Figure H shows the results for [i = — 1.8t, i.e., n » 0.41. 
As shown in our previous paper ||, the critical field 
increases markedly near fi = —2t due to an effect of 
the structure of the Fermi-surfacc in the present model. 
Therefore, for /i = — 1.8t, the critical field is very large 
at T — 0. It is found that the enhancement occurs espe- 
cially at low temperatures, where ifq — is the optimum 
direction. At rather higher temperatures, <p q = 7r/4 is 
the optimum direction, which gives critical fields slightly 
higher than that of </? q = 0. As shown in Fig. |[ other 
directions are optimum near the kink, where the curves 
of (fq = and ir/4 cross, although the differences in 
the magnitudes of the critical fields are very small there. 
Such behavior does not occur for fi = — 'it and // — > —At. 

Figure ^| shows the critical field for fi = — 0.5t, i.e., 
n w 0.77, where the Fermi-surface has a nearly square 
shape. It is found that the optimum q is oriented to 
the direction of the crystal axis (<^ q — 0) in a whole 
temperature region, and thus the kink does not occur. 
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FIG. 3. The solid and dotted lines show the critical fields 
for tfq = and 7r/4, respectively, when fi = — 1.8t, i.e., 
n ~ 0.41. In the inset, the solid and dotted lines show the 
temperature dependences of q = |q| along the transition line 
for ifiq — and 7r/4, respectively. 
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FIG. 4. Behaviors of the critical fields and q — |q| near 
the kink, where the lines of <^ q = and 7r/4 cross, when 
H — — 1.8t. The dotted and dashed lines show the critical 
fields of ifiq — and n/4, respectively. The solid line shows 
the final result of the critical field obtained by optimizing both 
the direction and the magnitude of q. The dot-dashed line 
shows the optimum ifiq. 
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FIG. 5. The solid, dotted, short dashed, dashed and 
dot-dashed lines show the critical fields for <^ q = 0, n/20, 
7r/10, 37r/20, and 7r/4, respectively, when jj, = — 0.5t, i.e., 
n ~ 0.77. The inset shows the temperature dependence of 
q = |q| along the transition line. 

In conclusion, we have studied the upper critical field of 
the superconductor on the square lattice in parallel mag- 
netic fields. In particular, we have examined anisotropy 
effects of Fermi-surface on the temperature dependence 
of the upper critical field. By changing the hole con- 
centration, we have controlled the shape of the Fermi- 
surface. The results are summarized as follows: (1) The 
kink of the critical field as a function of the temperature 
disappears due to an anisotropy of the Fermi-surface of 
some extent (Fig. ||); (2) When the anisotropy is weak, 
the kink remains (Fig. 0); (3) Depending on the shape 
of the Fermi-surface and the temperature, the optimum 
q is oriented to other directions than the symmetric di- 
rections ifiq = and 7r/4 (Fig. ||); (4) The critical field 
increases rapidly especially at low temperatures (Fig. H) 
near /i w — 2t, where the nesting is most effective ||. 
From the result (1), the absence of the kink in the exper- 
imental data which exhibits the upturn at low tempera- 
tures does not exclude the possibility of the FFLO state 
of anisotropic singlet pairing. 

When the orbital pair-breaking effect remains to some 
extent, for example, for large inter-layer coupling, the 
present result must be taken in another way. In this 
case, the direction of q coincides with the direction of 
the vortex lines for the coexistence of the vortex state 
and the FFLO state. Hence, ifiq must be regarded as the 
direction of the applied magnetic field in Figures ||, |^ and 
0. The temperature dependence of the upper critical field 
changes qualitatively as shown in those figures, by change 
of the direction of the applied field between (1, 0, 0) and 
(1, 1,0) directions, when the vortex state coexists. 
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